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Abstract 

We compute the complete one-loop finite temperature effective potential for elec- 
troweak symmetry breaking in the Standard Model with a Higgs potential supple- 
mented by higher dimensional operators as generated for instance in composite Higgs 
and Little Higgs models. We detail the resolution of several issues that arise, such as 
the cancellation of infrared divergences at higher order and imaginary contributions to 
the potential. We follow the dynamics of the phase transition, including the nucleation 
of bubbles and the effects of supercooling. We characterize the region of parameter 
space consistent with a strong first-order phase transition which may be relevant to 
electroweak baryogenesis. Finally, we investigate the prospects of present and future 
gravity wave detectors to see the effects of a strong first-order electroweak phase 
transition. 



1 Introduction 



The baryon asymmetry of the universe remains a mystery. Many ideas have been formulated 
in the literature, yet much uncertainty remains as to how the baryon asymmetry could arise. 
It is not even clear at what scale the initial asymmetry is produced. The Sakharov conditions 
for baryogenesis are baryon number violation, C and CP violation, and a departure from 
equilibrium. The Standard Model (SM) does not exhibit these conditions at nearly the 
strength required to produce the observed asymmetry given our standard cosmological 
assumptions, and thus it is expected that we must go beyond the SM in order to explain the 
asymmetry. 

The last Sakharov condition, departure from equilibrium, implies the necessity of a 
strong first-order phase transition. Since we know that the electroweak symmetry must 
be broken it is tempting to assume that the corresponding phase transition can satisfy this 
condition. As noted above, the SM is inadequate, but how far beyond the SM must one 
go to find the necessary out of equilibrium dynamics? This question has been addressed 
by a number of authors (e.g., see pQ for studies of the dynamics of the electroweak phase 
transition in various recent models). In Ref. [2J it was shown that if the Higgs potential 
is augmented merely by a H 6 operator, it can generate a strong first-order electroweak 
phase transition. As one can intuit, the scale suppressing this non-renormalizable operator 
must be in the neighborhood of the electroweak scale in order to generate a substantive 
effect on the phase transition dynamics. A tree-level analysis of this theory was conducted 
in [2J, with some further refinements in [3J, and it was concluded that a strong first order 
phase transition is possible even with a Higgs boson as massive as 200 GeV. Of course, for 
the presence of this H 6 operator to be compatible with electroweak (EW) precision data, 
a higher scale should suppress other dimension six operators, in particular those leading 
to oblique corrections. The analysis of [I] shows that the low energy effective theories of 
strongly interacting models, where a light composite Higgs emerges as a pseudo-Goldstone 
boson, have precisely this structure and single out H 6 as one of the dominant dimension 
six operators Q being suppressed by the decay constant of the strong sector, parametrically 
lighter than the cutoff scale of the model. For a fixed value of the strong decay constant, the 
compatibility with precision E W data is ensured by pushing the masses of vector resonances 
above 2.5 TeV [4j. In that case a decay constant as low as 300 GeV would be compatible 
with precision measurements. Our analysis should also apply to study the dynamics of 
electroweak symmetry breaking in Little Higgs models, and to the more general cases where 
the if 6 operator is generated by integrating out a heavy massive scalar field. However, some 
extra fine-tunings might be needed in that case to evade EW precision data. 

In this publication we extend the results of [21 [3J in several ways. First, we re-analyze 
the theory using the full finite temperature effective (nonrenormalizable) Higgs potential at 
one-loop. Second, we study the nucleation of broken phase bubbles and consider the effects 

-'in the case of a strongly interacting light Higgs boson, the general effective lagrangian includes four 
operators that are genuinely sensitive to the strong dynamics jl], i.e. suppressed by the strong decay constant 
and not the masses of the heavy resonances or the cutoff scale of the strong sector. In this context we can 
concentrate on the H 6 operator since it is the only one that affects the shape of the potential at tree-level, 
and thus it has significant effects on the dynamics of the phase transition as we shall see. 
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of supercooling on the electroweak phase transition (EWPT) within this more complete 
analysis. This is an important dynamical consideration of the phase transition that can in 
principle have dramatic consequences to when (and if) the phase transition happens. Finally, 
we investigate whether or not the gravitational waves emitted at the nucleation time can 
be detected by present and future interferometry experiments, which would provide another 
way to study the origin of EW symmetry breaking and another way to test the composite 
nature of the Higgs. We consider each of these points in the following three sections, and 
then make some concluding remarks. 

2 One-Loop Finite Temperature Effective Potential 

Once non-renormalizable interactions are allowed in the theory, as in our case, complete 
renormalization requires that the infinite set of higher-order operators be considered. How- 
ever, one is able to truncate the list of needed operators in a perturbative expansion of 
the inverse cutoff scale. To study the effect of new physics on the Higgs potential in this 
effective field theory context, it is sufficient to work at the order A -2 where A is the cutoff 
scale suppressing the effective operators. Higher dimensional operators will be sufficiently 
irrelevant to our problem and can be ignored. 

Our analysis is focussed on operators that affect the Higgs self-interactions. These effec- 
tive interactions parametrize the new physics responsible for EW symmetry breaking that 
become fully dynamical at about the scale A. Thus they can be used to generically constrain 
beyond-the-SM physics affecting the Higgs sector. Though EW precision measurements put 
severe constraints on the set of operators affecting the weak bosons' polarization tensors, the 
effective Higgs self-interactions are almost completely free parameters since the Higgs sector 
has not yet been probed directly by experiment. Thus the scale suppressing the operator 
if 6 we will focus on can be significantly lower than the cutoff scale of the (strongly coupled) 
model. This is in particular the case of composite Higgs models when the Higgs emerges from 
a strongly-interacting sector as a light pseudo-Goldstone boson 0j. The scale suppressing 
the H 6 operator is then /, the decay constant of the strong sector, a quantity 4-7T smaller 
that the cutoff scale. 

We start with the following classical effective potential for the SM Higgs [5]: 

V{H) = m 2 \H\ 2 + A|#| 4 + k\H\ 6 (1) 

where H T = (xi +%,V 9 + ^X3) /v2 which develops a vacuum expectation value (VEV) 
equal to Vq ~ 246 GeV. n~ l l 2 is identified at tree level with the decay constant of the strong 
sector - the details of this identification at one-loop are described later. We choose a vacuum 
configuration where only the real part of the neutral component has a constant background 
value: (p = <fi + h. The physical Higgs boson is h, and we use the traditional background 
field method [6J to evaluate the quantum potential for <p at one-loop. We focus on the main 
relevant contributions coming from the SU{2) L x U(l)y gauge bosons, the top quark, and 
the Higgs and Goldstone scalars. 
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As we briefly review in Appendix [Cj the quantum potential for the background value up 
to one-loop order at finite temperature in the Landau gauge (where ghosts decouple) is 



V eff {<f> t T) = V tree (4>) + AVi(4>, T) (2) 

with 

Vtreefr) = + \<f + fa, (3) 

AV^T) = E /tC^^ + ^ + ^W] (4) 

i=h,x,W,Z,t n=-oo J \ ) ~ 

where Ize — (co n , k) is the euclidean loop 4-momentum, uj n are the Matsubara frequencies in 
the imaginary time formalism, where u n = 2nirT for bosons (periodic on the euclidean time 
circle) and uo n = (2n + l)irT for fermions (anti-periodic on the euclidean time circle). The 
numbers of degrees of freedom for the relevant fields are n{h, x ,w,z,t} = {1,3,6,3, —12}. We 
include the fermion-loop minus sign in the definition of n t . 

Note that in the Landau gauge one must count all three degrees of freedom of each 
massive vector boson and the one degree of freedom of each Goldstone scalar. This may be 
qualitatively understood be recalling that the Xi Goldstone fields are independent quantum 
fluctuations away from the zero-temperature minimum. We present a quantitative argument 
showing this is not double counting in Appendix [Cj 

We obtain the background-dependent masses appearing in ^ by expanding the theory 
about the background value and reading off the quadratic terms for the various quantum 
fluctuations. In our dimension-six model the masses are 

m 2(0) = m 2 + 3X(j) 2 + ^K<f)\ (5) 

m 2 x ((f>) = m 2 + A0 2 + ^0 4 , (6) 

rn 2 w {j>) = mM = = f <A (7) 

where g,g' and y t are the SU(2) L , U(l)y and top Yukawa couplings respectively. At the zero- 
temperature minimum one recovers m\(vQ) = ml and m 2 (fo) = 0. Note that the expressions 
for the masses of the weak bosons (from the Higgs kinetic term) and the top quark (from 
the Yukawa coupling) are unchanged compared to the SM, and ([7]) are written to confirm 
our conventions. 

The one-loop correction ^ splits into a zero-temperature part and a T-dependent part [TJ 
[S] which vanishes as T — >• 0: 

AV^T) ee AVfit) +AV 1 T (0,T) (8) 
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with 



Atf(0,T) 



E 



n,- f d A ki 



%=h,x,W,Z,t 



(27T)' 



log [4 + 



2tt 2 



4 /"oo 



cWcfc 2 log 



(9) 
(10) 



AV\°(</)) is precisely the ordinary zero temperature effective potential, as it must be to be 
consistent since AV^((f),T) — > as T — > 0. The T = part, being UV-divergent, will be 
considered first in order to properly determine the renormalized parameters of the quantum 
theory. The finite temperature corrections will be treated afterwards. 



2.1 Zero Temperature Corrections 

At zero temperature the correction ^ reduces to the first term of ([8j, 



AV 1 (<j>,T = 0)=AV 1 °(<j>) = 



64tt 2 

i=h,x,W,Z,t 



uv 



which has been regularized in 4 — e dimensions, Cj = 5/6 (3/2) for gauge bosons (scalars 
and fermions) and Cjjv = ~ — 7s + log47r + 0(e). 

We work in the MS scheme to renormalize and evaluate our potential (see the Appendix [X| 
for an alternative, but ultimately equivalent, on-shell scheme approach). The full one-loop 
effective potential is 

Veff(<f>) = ^ + + + A^V) (12) 

where the parameters of this potential (m 2 , A, k) are bare parameters, but an implicit 5Vct 
will cancel their infinite pieces, leaving the finite pieces as the renormalized parameters. 

To determine the parameters of the lagrangian in terms of physical quantities, we must 
impose renormalization conditions at some chosen scale /i*. The renormalization conditions 
are 

VZ ff (<f> = vo,n*) = (13) 
V: ff (<P = v ,^) = m 2 h (14) 
^7/(0 = ^0,^) = e (15) 

The left side of each equation is the theory computation, and depends on the parameters 
of the theory (m 2 , A, k). The right side of each equation is a measurement (m^ and £) or 
related to a measurement {V'(vq) = is a requirement that the potential is at a minimum 
which recovers the correct Z boson mass). The VEV depends on the choice of scale as well. 
We define v$ to be equal to the VEV of the Higgs field in the Landau gauge at \i = mz such 
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that the MS Z mass is recovered. Performing our computations with the latest electroweak 
precision measurements j9], we find Vq = 246.8 GeV to a good approximation for a Higgs 
mass in our range of interest (115 GeV < < 300 GeV). This Higgs VEV is close to the 
246.2 GeV value in HO]. 



We can invert these equations to obtain the theory parameters as a function of measure- 
ments: 

ml = m 2 (m 2 h ,^v ,^) (16) 
A* = \(m 2 h ,£,v ,[M # ) (17) 
ft* = K(m 2 h ,£,v ,iJ,*) (18) 

Note, the parameters have scale dependence, and we have defined ml = m 2 (/i*), etc. 

Up to now we have glossed over some important subtleties. The physical Higgs mass must 
be defined at p 2 = ra|, whereas the one-loop effective potential is constructed for p = 0. To 
take account of this, and retain the label m\ for the physical Higgs boson mass, we need to 
rewrite the renormalization condition as 

m£-m£-£(m£) + £(0), (19) 

where £Qo 2 ) is the two-point function of the Higgs boson (numerically, we used the Loop Tools 
software [UJ to evaluate this two-point function). This approach has the added benefit that 
the IR singularity in V"ff(vo) as the Goldstone mass goes to zero is canceled by the IR 
singularity in E(0). We discuss these IR singularity issues in more detail in the Appendix |Bj 



The physical parameter £ is not a unique choice for how to parametrize the measured 
tri-Higgs coupling, and we wish to rewrite it in a more convenient manner. First, like the 
Higgs mass, the Higgs tri-scalar coupling has IR divergences at p = when the Goldstone 
bosons become massless. These IR divergences are also not dangerous because they are 
matched by the IR divergences of V"L(v ), and cancel in measured cross-sections. Thus, it 
is convenient to separate out this IR divergence when parametrizing the tri-Higgs coupling 
observable: £ = £p + ^m, where T IR contains IR sensitive Goldstone term^J Furthermore, 
since the tri-Higgs coupling £ in the SM is fixed once the Higgs mass is known, we would 
like our convention to reflect this manifestly in the decoupling limit of k — > 0, 

limf ^f SM = f| AJ + rf* f (20) 
For finite values of k, the deviations of from £j, M can be defined by convention to be 

7 



U^t S F M + 6 4 (21) 



This convention (i.e., the factor of 6) ensures that k x I 2 can be identified directly as the 
decay constant of the strong sector, /, at tree level. Putting these elements together, we can 



Explicitly, T IR is given by T IR 



:S2tt 2 



3ml(v )" m l(v o yiogml(v ) + 



where mf,($>) = ml + 



+ |k*</> . I/jj is given by the above expression where the limit — > is taken in mt.((j)). 
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now rewrite the third renormalization condition as 



vg f (vo) = z = el A/ + + T IR . (22) 



We emphasize that eq. (22) is merely a reparametrization of the tri-Higgs physical observable 
in terms of the decay constant, /, rather than £ for the benefits described above, and that 
£f M is a computable function of m h . 

Following the prescription provided above, all the parameters of our Higgs potential 
(m 2 , A, k) can now be written in terms of physical observables (vq, m^, f). Thus, we are now 
able to analyze the potential using physical observables as inputs. 



2.2 Finite temperature corrections 



From the splitting (|8j) of the full one-loop effective potential into a T = part and aT^O 
part, we get that the latter finite temperature component is: 



rp4 roo 

Aif (0,r) = E / dkk2 ^ 



i=h, X ,W,Z,t 



(23) 



" 2^ 2 vr 2 Jb I T 2 J + 2^ 2 tt 2 J/ I T 2 

i=bosons j=fermions 

where the upper (lower) sign stands for bosons (fermions). In the high-temperature regime 
(T ^> mj(0)), the Jj function expansions are 

7T 2 7T , /0 X 2 , X „„ / , , X 3 / 2 



■M*) = ^-V /2 -^log- + OU 3 log^— (24) 
z^o 12 o 32 at V est. 



7T 2 X 2 , X >_ f q X 3 / 2 



'/<*> =. -^-i log ^ + ( l3log ^rl (25) 



with logafe ~ 5.4076 and loga^ ~ 2.6350. Note that in [2] only the first terms in (24) and 



(25) were retained, which leads to the following approximate thermal one-loop correction: 



ATr <r . n-i \ sr^ niT 2 mf((b) v-^ nfT 2 m 2 (d>) 1 ~ , 9 , . 

AV! T G5 ^(0^)= E 24 + ^ 48 ~2 cT(j) +•••' ^ 



j=bosons j=fermions 

with c = (4m 2 /*; 2 + 3g 2 + g' 2 + 4y 2 - 12v 2 /f 2 )/16. 



The dominant contributions gathered in ( |26[ ) are simply a (positive) thermal mass which 
(meta)-stabilizes the origin of the potential at high temperature. This approximation was 
sufficient in |2j, and further refined in [3], to demonstrate the possibility of a strong first 
order PT within an effective extension of the SM. Fig. shows the discrepancy between 
the complete thermal correction and the high-temperature expansion around the critical 



temperature, illustrating the worthwhileness of using the integrals of (23) for the more 
detailed analysis. 
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Figure 1: Different potentials close to the critical temperature for rrih = 115 GeV and 
/ = 620 GeV (/ is the decay constant of the strong sector the Higgs emerges from). The 
dashed curve is the potential of [2] which includes only the thermal mass term of the Higgs, 
while the solid and dotted ones represent the full one-loop potential with (solid) and without 
(dotted) the ring diagram contributions. In blue, we have also plotted the imaginary part of 
the full one-loop potential with the ring contributions (solid blue) as well as the imaginary 
part of the ring contributions alone (dashed blue). This illustrates the cancelation of the large 
imaginary parts between the ring and the one-loop contributions, while there still exists an 
additional and smaller imaginary part for some values of due to a negative quartic coupling 



(see the discussion in Section 2.3.2 for details). An imaginary part of the potential can be 
interpreted as a decay rate of some quantum states of the scalar fields to some others but 
the imaginary part of the full potential is always tiny compared to the real part around the 
transition temperature and the system is stable enough throughout the entire time of the 
transition. 



2.2.1 Breakdown of perturbation theory and ring diagrams 

In thermal quantum field theory, the traditional perturbative expansion in terms of small 
coupling constants breaks down due to IR-divergences (inherent in massless models) gener- 
ated by long-range fluctuations appearing as soon as one moves to finite temperature [T2"] . 
For instance, taking massless A0 4 theory at finite temperature, one can show that the self- 
energy, which goes like A at first order, receives a subleading A 3//2 correction and not A 2 
as one would expect [13]. For our case, in the high-temperature expansion, or equivalently 



small mass expansion, of the thermal bosonic corrections (24), we also see a sign of this 
perturbation theory breakdown through the emergence of a monomial term of order 3/2. 
The main consequence is that, as it stands, we cannot trust the completeness of the one- 



loop result (23) because there are some higher-loop corrections of the same order [7], as if 
the effect of temperature is to "dilute" the one-loop correction to some multi-loop orders in 
the IR. Furthermore the leading part of these multi-loop corrections is all contained in the 
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Figure 2: Some generic examples of ring diagrams where each solid line may represent either 
a scalar, a fermion or a gauge field. The small loops correspond to thermal loops in the IR 
limit. They are all separately IR divergent, but their sum is IR finite. 



so-called ring (or daisy) diagrams shown in Fig. ^ . They are iV-loop diagrams where N — 1 
of them are "ring attached" to a main one. Since this "loop-dilution" is a finite temperature 
effect, the ring diagrams only need to be resummed in the IR-limit of vanishing momenta 
running in their petals [71] . It is also well-known that they can be taken into account by 
using propagators resummed in the IR [H] . By solving a Dyson-like equation, this turns out 
to simply shift the bosonic masses by a T-dependent constant as ml(<f)) — > m 2 {4>) + Ih(T), 
where Hb(T) is the self-energy of the (bosonic) field b in the IR limit, u — p — 0, known as 
a Debye mass (Ih(T) is labeled as lh(0) in [T3]). 

The higher-loop ring diagrams are needed due to IR divergences (i.e., m < T). On the 
other hand, the one-loop result is trustworthy for massive (i.e., m > T) particles, because 
the long-range fluctuations arising at finite temperature will never hit an IR mass-pole in 
such cases. Hence the ring diagrams will only contribute significantly at high-temperature 
(T/m — > oo) where the particles can be approximated as nearly massless. Also, this allows us 
to understand why only the bosonic degrees of freedom feel the breakdown of the perturbative 
expansion^] The reason is that only bosonic fields have a vanishing Matsubara frequency, 
recalling that u n equals 27mT for bosons and (2n + 1)ttT for fermions. Only this particular 
(zero-)mode will behave as a massless degree of freedom and generate IR-divergences at 
high-temperature, while the other (non zero-) modes u n act as a mass of order T and thus 
lead to negligible contributions. Therefore the fermionic propagators need not be resummed, 
because fermions do not have pole-mass in the IR. 

Applying the techniques of [H] to our theory, we compute the finite temperature mass 
shifts (Debye masses) that are needed in the ring diagram resummation: 

T 2 3T 2 v 2 

n h , x (T) = ^{m 2 h + 2m 2 w + m 2 z + 2m 2 )- — -± (27) 

Ti w (T) = 2 l<fT 2 (28) 

Hz(T) = f {ml - 2 m2w) T 2 -m 2 w {<P) (29) 

n 7 (T) = mU<t>) + f 4^ 2 - (30) 



3 In the gauge sector, only the longitudinal polarizations demonstrate this same breakdown of perturbation 
theory [H] . 
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Note that these II(T)'s are computed in the high-temperature limit of the unbroken phase 
which is justified by the ring diagrams being irrelevant for T < as we have discussed. 

At high temperature the photon and Z are not mass eigenstates, but one can treat them 
as mass eigenstates in this computation with the above-given Debye masses and obtain the 
correct resummed potential. 



2.2.2 Incorporating the ring corrections 



The traditional way the ring diagrams are implemented in the literature consists in shifting 
all the Matsubara modes for the bosonic fields. This is the so-called self-consistent method 
[15] where the potential Q is replaced by 



a Trself—con. 



>,T) 



E 



rijT 



+00 

E 



d 3 k 



log 



k 



+ ^ + m 2 (0) + n 4 (T) 



(31) 



The thermal shift of the gauge masses only for the longitudinal polarizations is understood, 
and n t (T) is simply zero. However, when applying this approach the UV divergent part 
becomes T-dependent through the n(T) and requires T-dependent counter-terms to be made 
finite. Indeed after doing to (31) the same splitting procedure we did to get (Tsl) , and after 



dimensionally regularizing the UV-divergent part, we get the following result: 



a TrO, self —con. 
" 1+ring 



E 

i=h, X ,W,Z,t 



K 2 (0) + n t (r)) 

64tt 2 



2 r 



log 



m 2 (0) + n,(T) 



— d — C 



uv 



(32) 



where the Cuv factor depends on T. This standard technique clashes with physical intuition 
since it would mean that the UV behavior of the theory depends on the IR dynamics. Al- 
though this mixing is not introducing any calculational errors to our working approximation, 
one can avoid it by simply shifting only the u n = Matsubara modes which carry the leading 
contribution from the ring diagrams relevant at one-loop order. 

As argued above, the dilution of the one-loop correction happens only for massless modes. 
Hence all the corrections we seek within the ring diagrams are gathered when resumming 
only the zero-mode of the propagator in the IR. Doing so, Q is to be replaced by 



AVi 



l-\-ring \}Pi 



T) 



E 

+ " 



E' 



d 3 k 



2 I ^ J (2vr) 3 

n=— 00 v ' 

log \k 2 + m 2 (0) + IU(T) 



log k 2 + ^ + m 2 (0) 



(2tt) 3 

= AT4(0,T) + Ay rinfl (0,T) 



(33) 
(34) 



where the prime means that the zero modes are excluded from the sum. We can easily 
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extract the ring part from the last expression and we find 



AV t 



ring yVi 



T) 



E 



i=h, x ,w,Za 



E 



i=h,x,W,Z,~i 



TijT 

Air 2 

127T 



dkk 2 log 



1 + 



n,(T) 



k 2 + m 2 ((f)) 



m^-Ko/o + n.cr)) 



3/2 



(35) 



where an irrelevant (infinite) constant has been ignored in the second line, and n{h, x ,w,z,j} = 
{1, 3, 2, 1, 1}. Notice that AV r i ng includes a monomial of order 3/2 which proves a posteriori 
the existence of a perturbation theory breakdown in evaluating the Higgs potential. Fur- 
thermore, these extra corrections modify the cubic term in m;(0), which partly controls^] the 
strength of the first order phase transition. Thus, the addition of these terms is critical for 
our analysis of the electroweak phase transition. 

In summary, the full T-dependent renormalized effective potential at one-loop is 



VeffW 



m 



* J.2 



+ 



A 



+ 



K 



* J$ 



i=h,x,W,Z,t 



i=bosons 



2tt 2 



J^2 



E 



64vr 2 
n,:T 4 



log^M-C 



j=fermions 



2tt 2 



-Jf 



J^2 



+ 



E 



n/T r 



i=h,x,W,Z,-y 



127T 



mf(0)-(m 8 2 (0) + n,(T)) 



3/2' 



(36) 



where definitions of all terms are given above. This is the potential we analyze for the 
remainder of the article. 



2.3 Reality of the quantum potential 

As the scalar masses become negative, the various contributions we obtained for the quantum 
potential develop some imaginary parts which we discuss below for both the T = and T^O 
cases. 



2.3.1 Imaginary part at T = 



In the zero-temperature limit, the logarithm of (36 ) leads to the following scheme- independent 
imaginary partj^] 

[AV?(<P)} = e(-^ 2 (0)) nil T: ( _ 0)|4 (37) 



i=h,x 



64tt 



4 The negative quartic coupling, of course, is another source of a potential barrier for the first-order phase 
transition. 

5 On the principal sheet, the imaginary part of the logarithm is taken to satisfy — ir < 3m log < ir. 
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where 0(— m 2 (</>)) is the Heaviside function which equals 1 when the field i is tachyonic, 
and zero otherwise. The Higgs boson can obtain a negative mass squared for some values 
of its VEV, originating from the fact that the classical potential is not convex everywhere. 
Indeed, depending on the cutoff value, either the origin is unstable (/ 2 > 3t>o/2m 2 ) or a 
potential barrier separates two local minima (f 2 < 3vQ/2m 2 l ), both of which lead to concave 
regions of the effective potential as a function of the VEV. A similar analysis shows that 
the Goldstone boson can become tachyonic for some values of the VEV as well, leading to 
another contribution to the imaginary part of the effective potential. However, we shall see 
shortly that the imaginary part (37) exactly cancels out with another contribution coming 
from the finite temperature corrections for the temperature range we are interested in for 
the phase transition. 



2.3.2 Imaginary part at T ^ 



At finite temperature both the integrals of (23) and the ring contributions (35) are spoiled by 



imaginary parts when scalar fields are tachyons. In the high-temperature limit, the imaginary 



part of (23) is (see (24)): 



3m 



T 



£e ( - 

i=h,x 



-m-(0))n. 



\mi{4>)\- 
64tt 



+ 



12tt 



(38) 



The first term cancels the imaginary part from the logarithm of the T = potential correction 
(37), while the second is only compensated when the ring diagrams are added, since their 



imaginary part is given by 



3m [AV ring ((f),T)] 



12tt' 



(39) 



as long as the temperature satisfies m 2 (0) + IT(T) > for all <fi. Although somewhat more 



complicated algebraically to show (see Appendix C.4 for details), this cancellation occurs 
also for smaller temperatures of order T ~ |mj(0)|. 

Nevertheless and despite this cancelation, the potential is not everywhere real because for 
some values of T and 2 , m 2 (0) + Ilj(T) < and the second term of the ring correction (35) 
becomes imaginary. In the SM this term does not lead to an imaginary part once the 
temperature (meta) stabilizes the origin since the SM scalars could only become tachyonic 
for a negative quadratic coupling in the Higgs potential. Thus, the SM potential is real 
as long as the origin is (meta)stable. On the other hand, with the additional if 6 piece in 
the potential, the scalar masses can be negative also through a negative quartic coupling, 
allowing this additional imaginary part to the potential at temperature around the critical 
temperature. 

An imaginary part of the potential can be interpreted as a decay rate of some quantum 
states of the scalar fields to some others [16J. Thus, one can rely on the real part of the 
potential as long as its imaginary part remains small enough to consider the field stable during 
the phase transition, in which case it can be discarded. We checked that the imaginary part 
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of the one-loop potential is always tiny compared to the real part around the transition 
temperature, thanks to the previously demonstrated cancelations of large imaginary pieces. 
Thus, we conclude that the system is stable enough throughout the entire time of the 
transition, and that its dynamics is driven by the real part of the one-loop potential we 
computed. 

3 Dynamics of the Electroweak Phase Transition 

Now that we have the formalism developed for our analysis of the finite temperature Higgs 
potential at one loop, we are in the position to study the dynamics of the phase transition. 
One of our first considerations must be the analysis of when (and if) the phase transition 
actually occurs. This is not simply a matter of determining the temperature at which the 
symmetry breaking minimum becomes the global minimum. An analysis of the energetics 
of bubble formation must be undertaken for a more complete picture. The nucleated 
bubbles can then undergo collisions and the surrounding plasma experience turbulence, 
which generate gravity waves that could possibly be detected in experiments. We discuss 
these issues in this section. 

Throughout this section, we report our numerical results of various relevant quantities 
as contour plots that scan the allowed region of the parameter space (m^,/). We recall 
that rrih is the physical Higgs mass while / is the decay constant of the strong sector (or 
more generally the energy scale suppressing the H 6 operator) physically defined through 
the triple Higgs self-interaction as defined in the previous sections, and we work in the MS 
scheme for \x = mz- The bounds delinating the region of first-order phase transition are 
both numerically computed using the complete one-loop potential at finite temperature. 
The lower one is set by requiring that EW symmetry is broken at T = and restores at 
high temperature, while above the upper bound the Higgs vacuum is likely to undergo a 
second-order phase transition or a smooth crossover. In general, determining the latter is 
not an easy task as it requires a non-perturbative analysis of the effective potential when 
the transition is not strongly first-order [T7j. Indeed, the phase transition always appears 
first-order at the perturbative level, even though very weakly. Moreover, as / increases one 
tends to recover the SM potential, which leads non-perturbatively to a continuous crossover, 
instead of a weak first-order transition at one-loop, for > 80 GeV [T8] . We estimated the 
upper bound by considering that as soon as the phase transition is as weak as in the SM for 
rrih = 80 GeV, it is likely to be a crossover. 

3.1 The onset of nucleation and EW baryogenesis 

The effective potential ensures the presence of a potential barrier at finite temperature which 
is a necessary ingredient to have a first-order phase transition. It proceeds by spontaneous 
nucleation of non-vanishing VEV bubbles into a surrounding symmetric metastable vacuum. 
As soon as the universe cools down to a critical temperature T c the symmetry-breaking 
vacuum becomes energetically favorable and then thermal fluctuations allow the bubbles 
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to form. However, the temperature of the transition is not necessarily close to T c . Once 
created, a bubble needs to consume a part of the latent heat liberated in order to maintain 
its interface with the symmetric phase surrounding it. It turns out that for T just below T c 
it is often the case that the bubbles are too small and surface tension makes them collapse 
and disappear. Hence the phase transition effectively starts at a smaller temperature when 
enough free energy is available to permit the nucleation of sufficiently large bubbles that can 
grow and convert the entire universe into the broken phase. This supercooling phenomenon 
can substantially delay the phase transition and thus modify the spectrum of gravity waves 
significantly, as we shall discuss shortly (important supercooling effects were also observed 
in some of the analyses of Ref . pQ ) 



3.1.1 When does the nucleation start? 

Although the probability to tunnel via the excitation of SU(2) instantons is very tiny, about 
exp(— 0(100)), the decay of the false vacuum can nonetheless proceed through thermal 
fluctuations which help to overcome the potential barrier. The rate per unit of space-time 
for this process is given in the semi-classical WKB approximation by T ~ c~ Se where Se 
is the euclidean action for the Higgs VEV evaluated on the so-called bounce solution of the 
euclidean equation of motion [19]. For temperatures much higher than their inverse radius, 
the bubbles overlap in euclidean time and feel the IR breaking of Lorentz symmetry [20, 21 J, 
in which case the bounce solution is 0(3)-symmetric and is the solution of 

d 2 <\> h | 2 #5 | dVj^T) _ Q 
dr 2 r dr dcfib 

subject to the boundary conditions 

fe ( r - oo) = and ^ = ° } = 0. (41) 

dr 

The bounce solution physically represents the Higgs VEV profile of a static unstable (either 
expanding or shrinking) bubble, and r measures the distance from the bubble center. For 
such a static solution of the equation of motion, the action factorizes as Se = S 3 /T, with 



Ss = / drAirr 2 



1 

2 V dr 



2 



+ V(<f> b ,T) 



(42) 



Moreover for small temperatures of the order of the bubble size, we replace the 0(3) bounce 
for the 0(4)-symmetric solution which minimizes the action when the breaking of Lorentz 
symmetry is not significant. Finally we use the traditional overshooting/undershooting 
method to numerically solve the equation of motion. 

There is a supercooling effect that can delay the onset of the first order phase transition 
to temperatures much smaller than 100 GeV. A first order phase transition can only proceed 
in the presence of a potential barrier separating the two vacua and the nucleation could 
potentially start at a temperature T n far below that of T c . This is especially likely in the 
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case where the barrier persists down to T = 0. Since the amount of supercooling is controlled 
by the size of the nucleated bubble, one needs to take into account that the phase transition 
proceeds in an expanding universe. One can thus consider that the nucleation starts at the 
time when the probability of creating at least one bubble per horizon volume is of order one. 
This condition guarantees the percolation of bubbles in the early universe and translates 
into the following criterion for determining the nucleation temperature: 

~ -4 log ( —) =s- ~ C(130 - 140) for T n ~ 100 GeV. (43) 

T n \mpi ) T n 

where mpi = Mpi/\f&K is the reduced Planck mass. 

The contours of constant nucleation temperature are reported in the left panel of Fig. [3] 
We point out that there exists a region (painted red in Fig. [3j with low / and m h < 225 GeV 



such that the criterion eq. (43) is not satisfied, meaning that the expansion of the universe 
does not permit the bubbles to percolate. Thus the nucleation never starts and the universe 
remains trapped in a symmetric vacuum. In addition, the right panel of Fig. [3] helps one to 
realize further the numerical significance of the supercooling effect by plotting the deviation 
of the nucleation temperature T n from T c . We see that, for large values of /, the deviation 
is not significant since the potential barrier disappears at a temperature not much less than 
the critical one. On the another hand, as soon as one lowers /, the barrier persists to lower 
and lower temperatures, making the supercooling delay of the phase transition important. 
Thus the knowledge of the nucleation temperature becomes necessary to clearly understand 
the dynamics of the phase transition in this region. 



3.1.2 Saving the baryon-asymmetry from wash-out 

Understanding the dynamics of the phase transition is a worthy endeavor on its own; however, 
one of the key reasons for understanding the nature of the EW phase transition is to 
determine if a baryon asymmetry can be produced and survive the process. Calculations in 
the previous sections enable us to refine some of the results of [2], where the possibility of a 
strong first order phase transition was first demonstrated. 

So far we have computed the crucial ratio (<f>(T)) /T at the nucleation temperature in the 
cases where only the thermal masses are included and where the complete one-loop potential 
is used. This allows us to compare the effect on the wash-out criterion of the supercooling 
of the phase transition and the usefulness of the one-loop potential. The contour plots of 
Fig. Q show the common fact that the lower the value of /, the stronger the phase transition 
for a fixed Higgs mass. The qualitative result of considering the temperature delay from T c 
to T = T„ is that for a given point in the parameters plane, the phase transition is generically 
stronger at T n . Indeed not only is the nucleation temperature potentially much smaller than 
T c , but also the value of the Higgs VEV grows as the universe cools down. 

Another important result for the baryon-asymmetry of the universe, is that it can be saved 
from the wash-out through sphaleron processes, namely {4>(T))/T > 1, for a not-so-small 
value of /. Indeed, in order to allow baryogenesis during the EWPT in the approximation of 
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Figure 3: The left panel of this figure shows contours of the nucleation temperature T n in 
the allowed region for an EW symmetry-breaking first order phase transition (/ is the decay 
constant of the strong sector the Higgs emerges from, and is the physical Higgs mass). 
Below the red lower bound the EW symmetry remains intact in the vacuum while above the 
blue upper one the phase transition is second order or not even occurs. Within the red band, 
the universe is trapped in a metastable vacuum since no expanding bubble is nucleated and 
the transition never proceeds. The contours are from left to right for T n = {50, 100, 150} GeV. 
The right panel of this figure shows contours of the relative deviation of the nucleation 
temperature from the critical one: 6t = (T c — T n )/T c . This measures the degree to which 
the phase transition is delayed by the overcooling effect. The contours are, from above, for 
e T = {10- 3 ,10- 2 ,0.1,0.3}. 

[2] some fine-tuning might be required in some approaches without any particular dynamics 
to make the suppression scale of the dimension six operator in the Higgs sector relatively 
smaller than the TeV scale required in the gauge sector to pass EW precision measurements. 
But the full one-loop potential tells us that for values of the Higgs mass above the current 
experimental bound / can be larger - as large as 1.2 TeV - and the baryon-asymmetry can 
still freeze out. 



3.2 Gravitational Waves 

As a bubble expands a part of the latent heat released accelerates the bubble wall and 
introduces turbulent motions in the hot plasma. After bubbles collide, spherical symmetry 
is broken. This enables gravitational radiation to be emitted. The turbulence of the plasma 
after bubble collisions is another important source of gravitational radiation (see [23] for an 
introduction to the physics of gravity waves). In the following, we characterize the spectrum 
of gravitational radiation that one can expect from the first order phase transition we have 
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Figure 4: Plot of the ratio £ n = (<f>(T n )) /T n characterizing the strength of the phase transition 
using the thermal mass approximation of [2] (left) and the complete one-loop potential 
(right). The contours are for £ n = {1,2,3,4} from top to bottom. / is the decay constant 
of the strong sector the Higgs emerges from, and is the physical Higgs mass. 



detailed in this article. We compare these results with the sensitivities of current gravity 
wave detectors, and of proposed gravity wave detectors of the future. 



3.2.1 Characterizing the spectrum 

Previous studies [24"l 1251 126] of the gravity wave spectrum culminate in showing that it can 
be fully characterized by the knowledge of only two parameters derived ultimately from the 
effective potential The first one is the rate of time-variation of the nucleation rate, named 
(3. Its inverse gives the duration of the phase transition, therefore defining the characteristic 
frequency of the spectrum. The second important parameter, a, measures the ratio of the 
latent heat to the energy density of the dominant kind, which is radiation at the epoch 
considered: a = e/p ra d. They are both numerically computed from the effective action S3/T 
at the nucleation temperature as follows. The time-dependence of the rate of nucleation is 
mainly concentrated in the effective action and (3 is defined by (3 = —dSE/dt\ . Using the 

6 This conclusion is valid under the assumption of detonation. However, in practice the bubble expand in 
a thermal bath and not in the vacuum and friction effects taking place in the plasma slow down the bubble 
velocity. Therefore, it might be important to consider the deflagration regime as in Ref. |27j . When the 
phase transition is weakly first order, we obtained under the approximations of |28j a wall velocity lower 
than the speed of sound. However, in the interesting region where the phase transition gets stronger, we 
approach the detonation regime and the approximations of [2Sj have to be refined to accurately compute the 
wall velocity. 
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Figure 5: The panel on the left contains contours of the latent heat a = 
{5.1CT 3 , 1CT 2 , 5.10~ 2 , 0.1, 0.5} from top to bottom. The panel on the right draws contours of 
the parameter, /3/H n , measuring the duration of the phase transition. From above one has 
(3/H n = {10 5 , 10 4 , 10 3 ,200}. / is the decay constant of the strong sector the Higgs emerges 
from, and mh is the physical Higgs mass. 



adiabaticity of the universe one obtain the following dimensionless parameter: 

d ( S 3 



dT \ T 



(44) 



where H n is the expansion rate when nucleation starts. The latent energy is the sum of the 
amount of energy AV seperating the metastable vacuum to the stable one and the entropy 
variation AS between these two phases. Hence one has: 



-AV - TAS 



-AV + T 



dV 
df 



(45) 



The left and right panels of Fig. [5] show contours of constant a and f3/H n , respectively, at 
the time of nucleation. 



3.2.2 Observability at interferometry experiments 

Future interferometry experiments could offer us a way to observe the EWPT. A detailed 
analysis of the potential to directly see gravitational waves from the first-order phase tran- 
sition can be compared with the sensitivity expected from the correlated third generation 
LIGO detector on earth and the LISA and BBO detectors in space. A general analysis that 
we utilize has been presented in [22], where both bubble collisions and turbulent motions 
were considered. Qualitatively, gravity-wave detectors will give us a better chance to observe 
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Figure 6: Example of gravity wave spectrum produced during the EW phase transition both 
by turbulence (left peak) and collision effects (right peak slightly emerging from the tail of 
the turbulence spectrum). This plot is for m h = 115 GeV and / ~ 600 GeV where a = 0.51, 
/3/H = 89 and T n = 39 GeV. Note that suitable values of a, /3/H to get a strong signal 
always imply a small nucleation temperature (< 100 GeV) due to important overcooling 
effects that drag the peak below the lower bound of the space-based detectors frequency 
band (~ 10~ 4 Hz), making the gravity waves delicate to observe. 

the phase transition today if the latent heat energy released is large and the emission lasts 
a long time. This can be understood easily by recalling that the power spectrum is given 
by the square of the quadrupole moment of the source which in turns scales as the kinetic 
energy over the time of emission In other words, typically a has to be 0(1) and /3/H 
as small as 0(100) to get a sufficiently high energy density flh 2 > 10 -10 . 

Relying on our effective (nonrenormalizable) potential approach, we find that generically 
the dynamics of the first order EWPT beyond the SM generate too weak gravity waves 
to observe except for a tiny region of the parameter space. Namely, by looking closely at 
Figs. [5] one can see that for a Higgs mass slightly above the LEP2 bound, > 115 GeV, 
and a relatively low scale, / ~ 650 GeV, we get at best a ~ 0.5 and (3/H ~ 100. The 
corresponding nucleation temperature in this region is about 50 GeV, according to Fig [3j 
For such a temperature scale, only LISA and BBO will be sensitive to the emitted spectrum 
of gravity waves, according to the results presented in Figs. 3 and 4 of [22]- Its detectability 
is probably beyond the capability of LISA. This result is in qualitative agreement with the 
results of [30]. Indeed LISA requires at least values of a > 0.6 for (3/H ~ 100 in order to 
see the characteristic peak from turbulence while the collision peak starts to be probed for 
a > 0.8. On the other hand, BBO should be able to observe both peaks if a is around 0.3 
(keeping (3/H ~ 100). 

Thus it seems that one will have to wait until the launching of the second generation 
of space-based interferometers to really study the EWPT through gravity wave detectors 
within this framework. Moreover this would be possible only in the maximizing case where 
the Higgs mass is close to its current experimental bound and the composite scale of the 
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Higgs is relatively low. 



4 Conclusions 

In this article we have reported on a complete computation of the one-loop finite temperature 
effective potential in models where the Higgs boson is composite and emerges as a light 
pseudo-Goldstone boson of a strongly interacting sector (our analysis could also be relevant 
for studying the dynamics of electroweak symmetry breaking in Little Higgs theories). These 
models are characterized by higher dimensional operators in the Higgs sector suppressed by 
the strong decay constant, /, a scale parametrically smaller than the cutoff of the strong 
sector. Interestingly, by following the details of the phase transition dynamics, the parameter 
space of a strong first-order phase transition has actually grown for large value of /, and 
shrunk for small value of / cutoff, compared to the tree- level result found in [2] . It has grown 
at the higher end by going beyond the high temperature approximation. The parameter space 
has shrunk on the lower end, since we found that bubbles cannot be nucleated well enough 
there to overcome the effects of an expanding universe. We encountered some subtleties 
along the way, including infrared singularities and imaginary components to the potential, 
that were resolved. 

It was also necessary to compute the details of the phase transition dynamics in order 
to investigate the possibility of detecting gravitational radiation from the first order phase 
transition occuring in the early universe. After bubbles are nucleated, their collisions and 
subsequent turbulence in the plasma give rise to gravity waves. In the assumption of a 
detonation regime, the effects depend on only two parameters, the latent heat a and the 
duration of the phase transition /3 _1 , both of which can be determined by solving the bounce 
equation, and analyzing the full one-loop finite temperature effective potential at the scale of 
the nucleation temperature. Although LIGO and LISA are likely not sensitive to these effects, 
we found that BBO, a planned second generation experiment of space-based interferometers, 
could be sensitive to the gravity waves produced during this phase transition. 
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A On-Shell Renormalization of the T = 



Potential 



The on-shell scheme identifies lagrangian parameters as physical parameter (i.e., observ- 
ables). It is the scheme employed by [3T], although we augment that discussion by describing 
a self-consistent approach with higher order operators, and describe the details of how IR 
divergences from massless Goldstone bosons cancel. 

Renormalizing our theory in the on-shell scheme is most convenient when we begin by 
writing the full potential in the following form: 
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(47) 
(48) 



which has been regularized in 4 — e dimensions, Cj = 5/6 (3/2) for gauge bosons (scalars and 
fermions) and Cjjy = - — 7^ + log47r + 0(e). In this parametrization of the tree-potential, 
the scalar (^dependent masses are: m\(<t>) = A(3</> 2 — Vq) + 3k(50 4 — 6t>Q0 2 + Uq)/4 and 



m 2 x {<j>) = A(0 2 - v 2 ) + 3^ 
expectation of the Higgs field, A 
below. 



,2\2 



/4. The on-shell scheme imposes that vq is the vacuum 
and k = l/f 2 . The precise meaning of f 2 is defined 
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The counter terms, <5q, are determined by the renormalization conditions: 
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where AS = E(m/ l )— E(0) and Ar = r(m/ l )— T(0) are needed to take us from the IR-sensitive 
and unphysical p = limit of the effective potential to p 2 = ml, where physical observables 
rrih and the tri-Higgs coupling £ p h ys are defined. Detailed computations demonstrating the 
cancelation of the IR divergences in this scheme are presented in Appendix |Bj 

We wish to have a more direct physical parameter that parametrizes deviations from the 
SM, and so we redefine 
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(52) 
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coupling in the SM is fixed with knowledge of rrih, and thus Cphy S * s determined completely 



which constitutes the definition of the physical observable /. Recall that the tri-Higgs 
coupling in the SM is fixed with knowledge 
by rrih and the other parameters of the SM: 
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(53) 



Since m 2 (vo) depends on 1/f 2 , this expression is technically equal to the SM one only in the 
limit of p —>■ oo, which is all that we need for the analysis to be self-consistent. 

We are now able to invert the renormalization conditions and compute the counter terms, 
which depend on the various derivatives of V®(4>), AE, and AI\ Upon expanding the result, 
one can express the renormalized full one-loop potential as 
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where all the /-dependence of the loop-order contribution to the potential is contained in 
the field-dependent masses, making the continuity of the decoupling limit explicit. 



B Cancelation of Goldstone Boson IR Divergences 

In this Appendix, we gather the detailed computations for the results mentioned in Ap- 
pendix [A] about smoothing the Goldstone IR singularity in the one-loop potential at zero 
temperature. First, we shall briefly recall how one moves from zero-momentum to on-shell 
scheme in the SM, as a warm-up for the dimension-six operator discussion that will come 
afterwards. 



B.l Review of the SM case 



In the SM the loop- integral of (47) can be renormalized by imposing the two conditions 
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which leads to the traditional form of the effective potential 
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where the scalar masses must be evaluated in the decoupling limit of the dimension-six 
operator (/ — ► oo). As we will review, m^o is an off-shell Higgs mass defined at p — 0. 
This is fine to use as long as no massless particle couples to the Higgs field |31j . If such 
particles like the Goldstone bosons are to be taken into account, one must move away from 
zero-momentum to avoid the pole-mass at p = that makes both m^ an d the one-loop 
potential IR divergent. 

In order to see how this can be done, we recall that near the symmetry breaking minimum 
(<f) = Vq) the renormalized effective potential can always be expanded in terms of lPI-Green 
functions evaluated at vanishing external momentum as follow: 



V eff (cf>, T = 0) = -J2 { ^^G^\p 2 = 0), (5f 



n=0 



where Q^ n '(p 2 ) are the n-legs renormalized 1PI Green functions for the physical Higgs scalar 
evaluated about the true vacuum (i.e., in the shifted theory). This expansion directly follows 
from the fact that the effective action may be intepreted as a generating functional of these 
1PI Green functions. Hence the second derivative of the effective potential at vq is simply 
the renormalized two-point function at zero-momentum: 



c?V e// (0,T = O) = _0( 2)(p2 = o)i (5Q) 
4>=vq 
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Given that the two-point function (the inverse propagator) of the Higgs is 

g {2 \p 2 )=p 2 -{ml R + nV 2 ))i (60) 

where m\ R and S(p 2 ) are the renormalized Higgs mass and one-loop Higgs self-energy, we 
see that imposing the renormalization condition (|56| leads to 



ml Q = ml R + Z(p 2 = 0). (61) 

justifying that mn,o is to be understood as the zero-momentum Higgs mass. In order to 
circumvent the IR divergences a natural choice would be to express the right-hand side of 



(56) in terms of physical parameters. The physical Higgs mass (mh) is defined as the pole 
of the one-loop resummed propagator (Q^ 2 '(p 2 = m\) = 0) and is given by solving the 
self-consistent equation: 

h,R 



m 2 h = m 2 +Z(p 2 = m 2 h ). (62) 



23 



P = 0\ 



-- h 



Figure 7: IR divergent diagram contributing to S(p 2 = 0). 



This allows us to rewrite ( 59 ) as 
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with AS = S(p 2 = m 2 h ) - S(p 2 = 0). 



When AS is absent the UV-finite one-loop correction at zero temperature for the SM 
has an IR divergent piece coming from the Goldstone contribution: 
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where m 2 = rrix SM \vo) will be kept non-zero as a regulator in what follows. Moving to the 



on-shell renormalization scheme (i.e., replacing (56) by (50)) results in the addition of the 
following term to V^f^: 



8Vi,{SM)(4>) 



AS (SM) 



2\2 



(65) 



Now T,( SM \p 2 = 0) receives an IR singularity from the diagram depicted in Fig. [?] which can 
be easily calculated to give: 
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Combining (|66|) with (|65|), one gets the following IR divergent contribution to the potential 
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(up to an irrelevant ^-independent term): 
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which exactly cancels out (64). 



B.2 Generalization to the non-renormalizable potential 

The presence of the dimension six interaction at tree-level forces us to set one more derivative 
of the potential to an extra measurable quantity. Focusing only on the decoupling limit this 
can be accomplished by 
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As can be easily checked, defining only the Higgs mass on-shell does not smooth out 
entirely the bad IR behavior of the one-loop potential in the non-renormalizable case. In fact 
/, being another parameter to be fixed at the quantum level, needs also to be renormalized 
away from zero-momentum to avoid the Goldstone pole, which is done by defining the 
renormalized three-point function of the Higgs boson on-shell. As for the two-point function, 



from (58) one gets 
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where pi denotes the external momenta of the three-point function. can be split into a 
tree-level coupling and a one-loop correction as 

£ (3) (ft 2 ) = -S3-r 3 (^) (7i) 

where #3 is the renormalized cubic self-couplings of the Higgs at tree-level. Similarly to the 



Higgs mass, we see that imposing (68) implies working with a parameter £o> or rather fo 



through (69), defined at zero- momentum which leads again to IR divergent behavior. We 
propose defining an on-shell cubic coupling at one-loop by[|] 



C phys = -g i3 \p: 



rrii 



93 + r 3 (p- 



m 



h)i 



(72) 



which translates into an on-shell (physical) definition of / by means of (52). Finally by 



of (51) 



expressing (68) in terms of physical parameters, we get the on-shell renormalization condition 

rf 3 K//(0,T = O) 



4>=vo 



93 + r 3 (p 2 = 

£,phys AIT^ 



0), 



(73) 
(74) 



where Ar = T^(pf 



mi 



-r 3 (rf = o). 



Now by enforcing the three renormalization conditions (49), (50) and (51) to set the 



counter-terms, we find that the zero-momentum potential is augmented by 



AT, 



o 2 ) + 



AS Ar 



16^ 



48ug 



2^3 



(0 2 --o 2 ) 



(75) 



We recall that in terms of m^o and £0 the effective potential develops a logarithmic IR 
singularity of the same form as in the SM but with /-dependent masses: 



AW) 



0(i\(IRdiv) 



n 



64tt 2 



x mHcfylogm 2 , 



(76) 



'Other physical definitions of the cubic coupling are possible, so long as they move away from zero- 
momentum to solve the IR issue. 
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Figure 8: IR divergent diagram contributing to T^{pf = 0). All the momenta are assumed 
to be zero in the external lines. 



while £o defined in (69) has a power-law divergence. Nonetheless T^{pj = 0) contains IR- 



divergent parts from the diagrams of Fig. M which are 



3n v mi ( _ 6v. 



r ( f Q) VR«») = gx_-^H h+J^ logmj (77) 



32n 2 m^Vy 



(7S 



Hence the power-law divergence of Ar cancels with the one of (69) making £, p h ys a well- 
defined quantity. The remaining logarithmic divergence of the three-point function along 
with the one from the self-energy, which turns out to be the same as in the SM, 

4 

n Tn 

Z {IRdm) (p 2 = 0) = \ogm% (79) 

gives after some simple algebra (up to irrelevant constant and 0(f~ 4 ) terms) 

5Vi{4>)V R ** = ^rritt) ^gm 2 x (80) 



which cancels with eq. (76). 



Finally, one finds that this procedure also leads to a UV and IR finite potential when a 
different higher derivative is chosen as a third renormalization condition. 



C Review of T ^ One-Loop Higgs Potential 

C.l The one-loop potential from the background field method 

The original method proposed by Jackiw in j6] to compute loop corrections to the classical 
potential is based upon expanding the action about (constant) background values for the 
various fields appearing in the theory. In our case only the neutral Higgs component has 
a non-vanishing VEV, and we recall here in a concise way how this method allows one to 
derive the one-loop correction given in (El). At the start, we consider the T = correction, 
and discuss the finite T corrections in the next subsection of this Appendix. 
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As an illustration we focus on a simple self-interacting scalar (real) field theory defined 
by the following generating functional: 



Z\j} = j[Vcl>]exp[i(S[<p]+j<p)] (81) 

where the notation ipj = j d A xip(x)j(x) will be assumed throughout this Appendix, and 
the action is S[ip] = f d 4 x[(d fl f) 2 /2 — Vo(y)]. Then one shifts the field by a x-independent 
background value (ip(x) = <fi + h(x)) where is assumed to be a classical field configuration 
and h represents a quantum fluctuation about it. We will now integrate out this fluctuation 
to get its effect on the potential up to one-loop order. To do so, one defines the shifted 
theory (whose dynamical field is now h) by expanding the action about its classical value: 



S[<l> + h]+j{<f> + h) = S{<j)] +j<j) + h 



6S 

5<p 



\ 1 5 2 S 

+ Ti^T—F— 

v=4> J 2 dip x d(f y 



<p=4> 



K + 



(82) 



where thanks to the equation of motion in the presence of a source the linear term vanishes. 
The ■ ■ ■ stand for higher (than quadratic) orders in h which lead to (at most) two-loop 
corrections [6]. One also easily obtains after an integration by parts that 



5 2 S 



tp=<f> 



-{U + V>i{4>))8\x-y). 



33) 



Plugging this expansion back into (81) one obtains 
Z[j] ~ jtsm+j® x f [ph]exp 



i, 5 2 S 
-h a 
2 2 



e i(sm+j<t>) x Det (□ + ^»(0))-3 . 



<p=4> 



34) 



We recall that by definition the effective action is the Legendre-transform of the logarithm 

S eff [<j>]=-ilogZ\j]-j<p, (85) 
which in our case, including the quantum fluctuations at one-loop, takes the form: 



S eff ^] = S[^] + -Tr\og 



5 2 S 



5ip x 8<p % 



tp=4> 



(86) 



Moreover S e tt can always admit a derivative expansion of the form: 



S, 



effl 



d 4 x [-K // (0) + A(0)(^) 2 + ---] 



^7) 



which defines precisely what one calls the effective potential. Since (f) is an homogeneous 
configuration in space-time, this simplifies to: 



Seff[4>] = -Weffffl, 
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where V is the volume of space-time we choose to keep finite for the moment. Besides, 
this homogeneity preserves the diagonality of (83) in momentum-space, which allows us to 
evaluate the trace in (86). This leads to the following expression for the effective potential: 

Veffffl = Votf) - '-V-'^ogi-k 2 + K'(0)) (89) 

k 

where the sum is over the eigenvalues of the □ operator in momentum-space. Finally by 
taking the limit of infinite space-time volume, one gets the well-known result: 

Veff(<f>) = V (<j>) l0 g(-^ + TO))" (90) 

The generalization for fields of higher spin that couple to <fi is 

K//(0) = W) + i V I tSj logdet (-iA(M)) , (91) 

j=flelds ^ ' 



where —iD(k,(j)) is the inverse propagator, 77 = —1/2 (1) for bosons (fermions) is the power 
of the functional determinant, and the det denotes an eventual determinant acting on either 
Lorentz or Dirac indices. 



C.2 Turning on the temperature in the effective potential 



The imaginary time formalism to go from quantum statistics at zero-temperature to thermal 
quantum statistics is by compactification of the euclidean time dimension on a circle of radius 
R = 1/2ttT. This correspondence is formally obtained in the path integral formulation of 
quantum mechanics |13j . However, it is worthwhile to give a quick intuitive argument. 



We begin with the generating functional for a scalar field in euclidean space-time (r = it): 



Z[j]= /[P0]exp 



(92) 



Now requiring the euclidean time to lie in the interval — 1/2T ^ r ^ 1/2T, and restricting 
the field to static configurations, one ends up with 



Z[j]= /[P0]exp 



1 

T 



(93) 



For vanishing source the space integral is nothing else but the energy (E[(j)]) stored in a 
(time-independent) field configuration 0, and the generating functional reduces to 



Z\] = 0] = / e 



E 



' ~ > e 

S'=all states 



-E S /T 



(94) 



which is the common partition function of statistical mechanics where describes all possible 
(static) configurations of a given system in equilibrium with a heat reservoir at temperature 
T. 
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Therefore the prescription to follow as soon as temperature is switched on is rather simple. 
It consists of Fourier expanding the fields among its eigen (Matsubara) frequencies u n and 
discretizing the imaginary time integrals by the following replacement rule: 



(95) 



For instance, applying (95 ) to momentum integral in (90 ) to implement the finite temperature 



correction leads to a potential of the form presented in Q: 



(27F logK + fc 



Km 



(96) 



C.3 Gauge degrees of freedom in Landau gauge 

Recalling the Goldstone equivalence theorem of gauge theory, one might doubt the necessity 
of counting the longitudinal polarization of a (massive) gauge field and its associated Gold- 
stone mode as independent degrees of freedom when computing the effective potential in the 
Landau gauge. Here we clarify this fact in the simple case of an abelian Higgs model. To do 
so, we explicitly compute the one-loop contributions of the [/(l)-gauge, ghost and Goldstone 
fields to the Higgs potential in the gauge. The effective potential turns out to be gauge- 
dependent, however there is no need to worry since it is not a physical observable. We work 
at T = but the following discussion can be driven the same way when the temperature is 
turned on, since we never evaluate momentum integrals. 

We begin with the gauge field (A^). It will affect the Higgs potential at one-loop through 
the following term: 

A^(0) = ~\ J logdet {-W-l{k)) , (97) 

where det acts on Lorentz indices. In the gauge the inverse propagator has the usual 
expression in momentum space: 

- iD~l{k) = (-k 2 + m\{<t>)) Xll v {k) + i (-k 2 + £m 2 A (0)) U^(k) (98) 

with n^(fc) = r]^ — k^kyjk 2 and U^ u (k) = k^k u /k 2 being the transverse and longitudinal 
projectors respectively. Since the traces of U T,L = 3, 1 and the determinant are invariants, 
we can move to a basis where the matrices II = CUC" 1 are diagonal and read: 

fl T = diag(0, 1, 1, 1) , fl L = diag(l, 0,0,0) (99) 

In this basis the determinant can be easily evaluated and gives: 

A^i A (0) = ~\ J 70i [Slog (-k 2 + m\{<t>)) + log (-k 2 + £m^(0)) + log£] (100) 
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Now we move to the Goldstone boson (x) an d ghost contributions which are 



AV? +aho *(<f>) 



d 4 k 



\og[-iD-\k))+i 



d 4 k 

Wf 4 



log (-iDj ost (k) 



Given that, in the abelian Higgs model, the inverse propagators are 



iD~\k) 



x 



A; 2 -mJ(0)-£m^(0), 



(101) 



(102) 
(103) 



with m x the mass the Goldstone receives from its Higgs couplings, we obtain 



AV{ 



X+ghost 



(0) 



d 4 k 



log (-k 2 + £m 2 A (4>)) + m - log ( 1 + 



m 2 J(j)) 



(104) 

Gathering (100) and (104) together, we see the first terms of each expression cancel out, 

m 2 {4>) 



leaving only (in euclidean space) 



AV 1 



A+x+ghost 



(0) 



d 4 kE 



31og(A;| + mi(0))+log 1 



k% + Zm 2 A {<t>) 



+ ■ 



(105) 

where • • • stand for constant terms irrelevant for the potential. Taking £ = to move to the 
Landau gauge, the last expression reduces to (up to an infinite constant) 



A yA+x+ 9 host^ 



d A k L 
W) 



[3 log (k\ + m 2 M) + log (k\ + m\ (0))] + 



(106) 



from which one clearly sees that, in this gauge, the factor of 3 for the massive gauge field is 
not altered by the addition of the Goldstone contribution. 

Another physically meaningful fixing choice is the unitary gauge £ — > oo. Sending 
the gauge fixing parameter to infinity in (105) implies the decoupling of the Goldstone 
contribution, as it should: 



AV 1 



A+x+ghost 



(0) 



d A k L 

W) 



[31og(4 + r4(0))] 



(107) 



Again the degrees of freedom of the gauge field are still 3 in this gauge. 



From this discussion we see that interpreting the factors in front of the log as the number 
of polarization states for the corresponding field is only (accidentally) true in both £ = 0, oo 
gauges. Indeed, if one takes for instance the 't Hooft gauge (£ = 1) the results are 



A+x+ghost 



(0) 



d 4 k L 

W) 



[2 log (k% + m? A {4>)) + log 01 + m 2 x {(j>) + m 2 A {(j>))] + 



where now the "degrees of freedom" of reduce from 3 to 2 by this gauge choice. 



(10* 
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C.4 Cancelation of imaginary parts at small temperature 



In Section 2.3.2 we have shown that in the high temperature limit, a cancelation occurs 
between imaginary parts of the one-loop potential and the ring corrections. Here we want 
to show that this cancellation occurs also for smaller temperatures of order T ~ |mj(0)|. 



Indeed by working out the integrals of (23) in order to isolate its imaginary part, we get 



4 l^iWI 

3m [Alf(0,T)] = ©(-"»i(0)xT / dxx 

i=h,x 71 J ° 

where n is a positive integer which ensures that 



(An + 1) 



7T 



J^2 



IT < 



71 



J"2 



+ 2n7T ^ 7T 



(109) 



so that one stays on the principal sheet when taking the imaginary part of the logarithm, 
whose branch is assumed to lie on the negative real axis of the complex plane. We can easily 
show that n = as long as T > T n = |m.j(0)|/37r, in which case the imaginary part becomes 



3m[A^(0,T)] = £ 6(-m?(0)^ /" * dxx 2 In 

i=h, x n Jo \ 



mi(4>)\' 

J"2 



x- 



6(-m 2 (0)K 

i=h,x 



|mi(0)| 4 + |mi(0)| 3 T 



64tt 



12tt 



(110) 



and reproduces the same cancellation with (39) and (37) as in the high temperature regime. 



For completeness we now consider the case of very low temperatures. As the temperature 
cools down below T n , one begins needing to shift the imaginary part of the log by multiples 



of 2tt to remain on the principal sheet of the complex plane. Furthermore, from (109) we 
see that 



2tt < (An + 1)tt 



\m(4>)\' 

Ji2 



x 2 < 2tt. 



Thus 



Q(-mmnJ^f^ < 3m [AV?(<f>,T)] < £ 6(-m 2 (0))^^#^ (112) 



6/T 



i=h,x i=h,X 

and we conclude that 3m [AV^((p, T)] vanishes as T goes to zero, as it should. 
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